Calculations of two-color interband optical injection and control of carrier population, 
spin, current, and spin current in bulk semiconductors. 
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Quantum interference between one- and two-photon absorption pathways allows coherent control 
of interband transitions in unbiased bulk semiconductors; carrier population, carrier spin polariza- 
tion, photocurrent injection, and spin current injection can all be controlled. We calculate injection 
spectra for these effects using a 14 x 14 k • p Hamiltonian including remote band effects for five bulk 
semiconductors of zinc-blende symmetry: InSb, GaSb, InP, GaAs, and ZnSe. Microscopic expres- 
sions for spin-current injection and spin control accounting for spin split bands are presented. We 
also present analytical expressions for the injection spectra derived in the parabolic-band approxi- 
mation and compare these with the calculation nonperturbative in k. 
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I. INTRODUCTION 

When a bulk semiconductor is simultaneously irradi- 
ated by an optical field and its phase-coherent second 
harmonic, quantum interference between one- and two- 
photon absorption pathways enables excitation of carrier 
distributions with interesting properties Such exci- 
tation, even without an external bias, can produce bal- 
listic photocurrents, 5 spin-polarized currents, 6 and pure 
spin currentsi^ Characteristically of quantum interfer- 
ence, these currents are sensitive to the phases of the 
two optical fields. In noncentrosymmetric semiconduc- 
tors, the phases can also be used to control the total 
population of photoexcited carriers^ and the net carrier 
spin polarization. 4,9 Which of these effects occur depends 
on the polarization states of the fields. 

These are examples of "n + m" coherent control 
schemes, in which a two-color light field controls a 
physical or chemical process by interference of n- 
and m-photon transitions^' 1 * 1 ' In semiconductors, 
"1+2" excitation has been discussed for impurity-band 
absorption^ free carrier absorption^" 1 ' quantum 
wellsj 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 and quantum wires'^ but our 
interest here is "1+2" coherent control of interband tran- 
sitions in unbiased bulk semiconductors^''' 2 ' 2 ' Such 
experiments have been performed with either (a) two 
fields, typically short pulses, one the generated second 
harmonic of the other, 2 i 4 i 5 i 6 i 7 i 8 i 9 i 27 i 28 . 29 i 30 i 31 i 32 i 33 i 34 i 35 or 
(b) a single ultrashort pulse having at least an octave 
bandwidth'-' 7 . 

Previous microscopic calculations of "1+2" processes 
in bulk semiconductors fall into two categories: ab initio 
density functional methods have been used for current 
injection 1 , and population control' while simple analyt- 
ical band models perturbative in k (with at most eight 
spherical, parabolic bands) have been used for current 
injectioni*2i 3 &224£ and spin-current injection' The for- 
mer are best suited for excess energies on the order of 
eVs, while the latter are only valid for excitation close 
to the band edge and cannot be applied to population 
and spin control, which vanish in such centrosymmetric 



models. 

In this article, we calculate "1+2" processes using an 
intermediate model that diagonalizes the k • p Hamilto- 
nian in a basis of 14 T-point states with remote band ef- 
fects included perturbatively. The model contains empir- 
ically determined parameters' 1 * 4 ' Fourteen-band mod- 
els (also called five-level models) have been used to cal- 
culate band structures ) 2i i 2' 4 '-' 4 ' linear" 4 ' and non- 
linear"'*"*' optical properties, and spin decoherence 
properties'" of GaAs and other semiconductors. Win- 
kler has recently reviewed 14-band models'*^ The model 
is nonperturbative in k and includes nonparabolicity, 
warping, spin-splitting, and interband spin-orbit cou- 
pling. We apply the 14-band model to the zinc-blende 
semiconductors InSb, GaSb, InP, GaAs, and ZnSe. 

We compare these results with analytic expressions de- 
rived in the parabolic-band approximation (PBA) based 
on an expansion in k about the T point of v nj7n (k) , which 
is the matrix element governing optical transitions. A 
one-photon transition is called "allowed" if the zeroth- 
order term in its expansion is nonzero, and called "for- 
bidden" otherwise. Two-photon transitions have two ve- 
locity matrix elements, and thus have a hyphenated label 
depending on the lowest-order terms in the expansions 
for each matrix element. For example, if both matrix 
elements are independent of k to lowest order, the two- 
photon transition is called "allowed-allowed" . For cur- 
rent injection and spin-current injection, we use expres- 
sions derived previously with an eight band model"i 
For population control and spin control, we derive ex- 
pressions with the 14-band model. 

The comparison between the PBA expressions and the 
numerical calculation establishes an important micro- 
scopic difference between current and spin-current con- 
trol on the one hand, and population and spin control on 
the other hand. Close to the band-gap, the former result 
from the interference of allowed one-photon transitions 
and allowed- forbidden two-photon transitions, whereas 
the latter result from the interference of allowed one- 
photon transitions and allowed-allowed two-photon tran- 
sitions. This difference was posited previously based on 
heuristic arguments "SI 
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Most of the early theory on semiconductor "1+2" 
processes processes conceptually separated the optical 
injection of densities and currents from the relaxation 
and transport of these quantities. We follow this ap- 
proach, and in this article, focus on microscopic calcu- 
lations of the optical injection. We note that relaxation 
and transport have been studied with an effective cir- 
cuit modelj^ii^i hydrodynamic equations^*^ Boltzmann 
transport in the relaxation time approximation^ a non- 
equilibrium Green function formalism^ and the semi- 
conductor Bloch equationsi 20 i 21 i 22 i 24 

We model the optical field as a superposition of 
monochromatic fields of frequency u> and 2uj: 

E(i) = E w exp(— iu)t) + E 2w exp(— i2ujt) + c.c. (1) 



and we sometimes write E,, 



E 



:/2u 



= \E„ 



/2u 



E 



and 



/2u 



exp(i(/) w /2w). We describe the fourteen- 
band model in Section [H] and use it to study "1+2" 
current injection in Section IIIII "1+2" spin-current in- 
jection in Section [V] "1+2" population control in Sec- 
tion ]W\ and "1+2" spin control in Section IVT1 We cal- 
culate the injection of each "1+2" process using micro- 
scopic expressions derived using velocity gauge (A • v) 
coupling in the long wavelength approximation, treating 
the field perturbatively in the Fermi's golden rule limit, 
and using the independent-particle approximationiSi^iS 
For spin-current injection and spin control, we use mi- 
croscopic expressions that include the coherence between 
spin-split bands. In Appendix 1X1 we justify the neglect 
of k-dependent spin-orbit coupling. The parabolic-band 
approximation results are derived and discussed in Ap- 
pendix [51 and compared with the numerical calculations 
in Sections IIIII I VII We summarize and conclude in Sec- 
tion Enj 



II. MODEL 

The fourteen-band model Hamiltonian, which in- 
cludes important remote-band effects to order fc 2 , and 
which we denote H14, is given explicitly by Pfeffer and 
ZawadskL^Si The fourteen bands (counting one for each 
spin) , are shown in Fig. ^ They comprise six valence 
bands (two each for split-off, heavy and light holes) and 
eight conduction bands (the two s-like ones at the band 
edge, and the six next lowest ones which are p-like). We 
now briefly review the derivation of -H14. 

The one-electron field-free Hamiltonian is H = Hq + 
Hso, where Hq — p 2 / (2m) + V, the potential V (r) has 
the symmetry of the crystal, and the spin-orbit interac- 
tion Hso is 



Elso 



h 



4m 2 c 2 



<r-(Wxp). 



where cr is the dimensionless spin operator, er = 
2S/H. Note that relativistic corrections proportional to 




r 6c — ri 



FIG. 1: A schematic diagram of the fourteen-band model, 
indicating band abbreviations (left), energies (center), sym- 
metry of the T-point states (right) , and one- and two-photon 
transitions. Te, IV, and Fg indicate irreducible representa- 
tions of the Td double group, whereas Fi and T4 indicate 
irreducible representations of the Td point group. Note that 
spin-splitting, nonparabolicity, and warping are not shown in 
this diagram. 



are Bloch states |nk) with energy hjj n {k). The asso- 
ciated spinor wave function 0„k (r) = (r|nk) can be 
written 4> n u_(r) — w„k (r) exp (zk • r), where the spinor 
functions u n k (r) have the periodicity of the crystal lat- 
tice. We use the notation |nk) to denote the kets 
for the u-functions; i.e. u„k (r) — (r|nk). Note that 
|nk) = cxp (— ik ■ r) |nk). The Hamiltonian for the it- 
function kets, known as the k • p Hamiltonian, is 57 - 58 



tt „ — ik r 

«k — e 



He lkr = H 



h 2 k 2 

2m 



fik- 



where the velocity operator v = (i/ti) [H, r] is 



1 

— P 



4m 2 c 2 



(cr x W) . 



(2) 



|er x Vy| have been neglected^ The eigenstates of H pointi 



The second term in v, the anomalous velocity, which 
leads to k-dependent spin-orbit coupling in H^, can be 
neglected for the processes we consider as shown in Ap- 
pendix ^ in the rest of this article, we assume that it 
vanishes. 

The states |n, k = 0) are a complete set of eigenstates 
for the Hamiltonian H on the space of cell-periodic func- 
tions. Thus cell-periodic eigenstates of can be ex- 
panded in the infinite set of states \n, k = 0) . The "bare" 
fourteen-band model truncates this expansion to a set 
of fourteen states, corresponding to the fourteen bands 
closest in energy to the fundamental band gap at the T 
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In a semiconductor of zinc-blende symmetry, the states 
{|n,k = 0) \n = 1..14} are conveniently expanded in the 
eigenstates of H , {\S} , \X) , \Y) , \Z) , \x) , \y) , \z)} <g> 
{IT); II)}) where, under the point group T^, \S) trans- 
forms like Ti, {\X) , \Y) , \Z)} and {\x) , \y) , \z)} trans- 
form like The {|T) > II)} comprises the usual spin 
1/2 states: 



(tk|T> = -Uk|4) = z 
(tkll) = ((lklt)r = x-iy. 

The non-zero matrix elements of (V7 x p) arc 

,4m 2 c 2 



(3a) 
(3b) 



{x\(yvx P y\z) = 

(x\ (W x p) v \z) = 

{x\ (w x P y \z) = 



3^ 
4rn 2 c 2 



A , 

An, 



4ra 2 c 2 . 

i A" 

3h 



cyclic permutations of these [e.g. (x\ {WV x p) y \z) = 
(z\ (W x p) x \y) = (y \ (W x p) z \x)], and those gen- 
erated by Hermitian conjugation of these. The above 
equations define the spin-orbit energies An and Aq, and 
the interband spin-orbit coupling A~i£2iSl The fourteen 
basis states {|n,k = 0) \n — 1..14} for P14 are 

\F 7v , ±1/2) = ±-L \Z) \a ± ) + \X ± *y) |a T ) (4a) 
|r 8 „, ±1/2) = \Z) \a±) + ±=\X± iY) jcfep) (4b) 
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(4c) 
(4d) 



|r 8 „,±3/2) = ±-^=LY±iY)|a±) 
|r 8c ,±l/2)=t|5)|a±> 
|r 7c , ±1/2) - ±-L |z) |a±) + |x ± *») |a T ) (4e) 

|r 8c , ±1/2) - T^f \z) |a±) + 4g I* ± iy) I«t> (4f) 



1 



|r 8c ,±3/2) =±^=|x±iy)|a±), 



(4g 



where = ||) and = \{). The states are la- 

beled with their transformation property under the dou- 
ble group for Td, and with a pseudo- angular momentum 
notation. In the basis (@J, Pk=o is diagonal except for 
terms proportional to A - . The connection between the 
eigenvalues of Pk=o for the T-point eigenstates and the 
eigenvalues of Pn is given by Pfeffer and ZawadskitM The 
nonzero matrix elements of momentum, which appear in 
Pk, are 

(S\ p x \X) = (S\ p y \Y) = (S\ p z \Z)= imP Q /H (5a) 
(S\ p x \x) = (S\ P y\y) = (S\ p z \z)^ imP'jn (5b) 

(x\ p y \ Z ) = (y\ p * \x) = (z\ P * \ y ) = (z\ p y \x) 

= (Y\p x \z) = {X\p z \y)^imQ/h. [ °> 



Eq. (0 defines the parameters Pn, Pg, and Q. They are 
sometimes expressed as energies Pp, Pp/, and Eq with 
the connections Pp = 2toPq/S 2 , etc. 

The "bare" fourteen-band model has eight empirical 
parameters E g , A , E' Q , A' , A - , P , Q, and Pq. Its 
quantitative accuracy is improved by adding important 
remote band effects to order k 2 using Lowdin pertur- 
bation theoryjSi which adds k-dependent terms to the 
truncated 14 x 14 Hamiltonian so that its solutions bet- 
ter approximate those of the full Hamiltonianj^i The re- 
mote band effects are governed by the parameters 71, 
72, 73, F, and Ck- The parameters 71, 72, and 73 are 
modified Luttinger parameters that account for remote 
band effects on the valence bands. They are related to 
the usual Luttinger parameters jil, 72L7 and 73L by the 
couplings with Tg c , T-? c , and r 8c bands, which are already 
accounted for in the "bare" fourteen-band model—: 



7i = 7ii 



72 = 72L 



73 = 73L 



Pp 


E Q 


3P S 


3P£ 


Pp 




6E g 


6E> 


E P 


Eq 


6E g 


6P 



E c 



3E' 



A' 



The parameter P accounts for remote band effects on 
the lowest conduction band, essentially fixing its effec- 
tive mass to the experimentally observed value. Finally, 
the parameter Ck is the small k-linear term in the valence 
bands ML The remote band effects can be removed by set- 
ting 71 = —1 and 72 = 73 = P = Ck = 0. The model 
includes neither remote band effects on the uc bands, 
nor remote band effects on the r6 C -r 8 „ and Tec-I^,, mo- 
mentum matrix elements, although such terms exist in 
principled 

In summary, P14 is a fourteen-band approximation to 
Pk that incorporates some remote band effects. It can be 
found in Eq. (5) of Pfeffer and Zawadzki, although with 
a slightly different notation. 4 -i With their notation on the 



left, and ours on the right: Pq 



-Eg, Pi 



E' n — E„ 



Ai = Aq, A = A", Pi = P^. Also, our A differs 
from theirs by a minus sign. Other authors have also 
used different notations.- 2 The fourteen bands are shown 
schematically in Fig.^along with the symmetry notation 
of the T-pomt states, and the notation used to label the 
bands. 



A. Material parameters 

Numerical values for the thirteen parameters of the 
model are listed in Table [Q for InSb, GaSb, InP, GaAs, 
and ZnSe. They are taken from the literature, where they 
were chosen to fit low-temperature experimental data. Of 
the two parameter sets discussed by Pfeffer and Zawadzki 
for GaAs, we use the one corresponding to a = 0.085 
that they find yields better results^i For InP, GaSb, and 
InSb, we use parameters from Cardona, Christensen and 
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TABLE I: Model parameters. 





GaAs 


InP 


GaSb 


InSb 


ZnSe 


E g (eV) 


1.519 


1.424 


0.813 


0.235 


2.820 


A (eV) 


0.341 


0.108 


0.75 


0.803 


0.403 


E' (eV) 


4.488 


4.6 


3.3 


3.39 


7.330 


A' (eV) 


0.171 


0.50 


0.33 


0.39 


0.090 


A~ (eV) 


-0.061 


0.22 


-0.28 


-0.244 


-0.238 


Po (eVA) 


10.30 


8.65 


9.50 


9.51 


10.628 


Q (eVA) 


7.70 


7.24 


8.12 


8.22 


9.845 


Pi (eVA) 


3.00 


4.30 


3.33 


3.17 


9.165 


JlL 


7.797 


5.05 


13.2 


40.1 


4.30 


72 L 


2.458 


1.6 


4.4 


18.1 


1.14 


■J3L 


3.299 


1.73 


5.7 


19.2 


1.84 


F 


-1.055 














C k (meVA) 


-3.4 


-14 


0.43 


-9.2 


-14 



FasaliSi For cubic ZnSe, we use the parameters given by 
Mayer and Rossler^ we use a calculated value of Ckr^ 
and we use A - = —0.238 eV to give a fc 3 conduction 
band spin-splitting that matches the ab initio calculation 
of Cardona, Christensen and FasaliSS Winkler used these 
same parameters for ZnSe, but took A - = 04^ There 
is more uncertainty in the parameters for ZnSe than in 
those for the other materials^ but we include it as an 
example of a semiconductor with a larger band gap. 

The parabolic-band approximation calculations use 
parameters from Table [J and average effective masses 
derived from the parameters in Table [I] 



Hamiltonian. But when remote band effects are included 
in a finite band model, they no longer hold. That is, 
v n m (k) calculated using © and eigenstates of H u is 
not equal to (nk\ Vk-Hi4 mk). We explicitly restore 
these identities by using (rtk| Vk-ffi4 |mk) to calculate 
Vn m (k). This approach can be described as including 
remote band effects in the velocity operator. It was used 
for an eight band calculation of linear absorption by En- 
ders et al^ This step is not critically important for the 
effects calculated here, since remote band effects are gen- 
erally small. 



C. k-space integration 

The optical calculations in this article have the form 
0, where 



9 = E / 6 ( k ) - 2huj ) > 



(10) 



where f cv depends on matrix elements and energies of 
eigenstates of H^, and where oj nm (k) = u> n (k) — ui rn (k). 
The integral in (fTUf) is understood to be restricted to the 
first Brillioun Zone, but we do not actively enforce the re- 
striction, since the photon energies considered here cause 
transitions well within the first Brillioun Zone. Writing 
k = {kevi $k, <Pk) m spherical coordinates, where k cv is 
the solution to 



B. Matrix elements 

The relations between matrix elements of the Bloch 
states and matrrix elements of the u-function kets are 

ffc 

v nm (k) = (nk\ v \mk) = (rik\ v |mk) H S nm , (6) 

m 

(nk| S |mk) = (nk| S |mk) (7) 

KUi 

(nkl v l S 3 Imk) = (nkl v l S 3 Imk) + (nkl S 3 Imk) . 

m 1 1 

(8) 

The matrix elements of the velocity operator, v, ne- 
glecting the anomalous velocity as discussed in Appendix 
lAl can be calculated using J2Jl, ©, and the right side of 
The matrix elements of the spin operator S, can be 
found from Eq. 0. The matrix elements of v l S 3 can be 
similarly found in the basis of eigenstates of Hq. Each of 
these can then be rotated to the basis (0} in which the 
states |mk) are expanded. 

It is well known that, in a crystal, v„„ (k) = Vk^n (k). 
More generally, 

v nm (k) = V k (nk| H |mk) = {nk\ V k ff k . (9) 

These identities can be proven from the definitions = 
e _lk ' r ffe !k r and v = (i/K) [H, r], even for a non-local 



f^Jcv (k C v, kl 4> k ) - 2hu = 0, 



(11) 



we have 



tt/2 r ir/2 



r-„J0 JO 



fc^sin^k/c (H u ) 



h(v cc (k) - v vv (k)) • k 



C?0kC?0k 



(12) 

where we have used Vw„ (k) = v, m (k) and the cubic 
symmetry of the crystal. It is numerically convenient to 
do the sum over any degenerate bands before the integral 
over #k and (/>k- 



D. Approximations 

The calculations of "1+2" effects in the following sec- 
tions are primarily labeled by the Hamiltonian used to 
approximate H^. The complete fourteen-band model is 
denoted H\^. The bare fourteen-band model, denoted 
^14-Baro, is -Hi4 without remote band effects. The 8x8 
subset of the fourteen band Hamiltonian within the basis 
{rg c , Tgu, r7„} is denoted Hs- The spherical eight-band 
model, denoted -f/sSph, is derived from Hg by setting 
C k = and replacing 72 and 73 by 7 = (272 + 373) 
it is a spherical approximation to the Kane model includ- 
ing remote band effects^ The aforementioned calcula- 
tions are non-perturbative in k; that is, in each case, the 
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Hamiltonian is solved numerically at each k. The per- 
turbative calculations of Appendix [B] are denoted PBA 
(parabolic-band approximation). 

The microscopic expression for each of the "1+2" ef- 
fects contains a sum over intermediate bands, which orig- 
inates from the two-photon amplitude. Unless other- 
wise noted, calculations include all possible intermediate 
bands (eg., Hn includes fourteen intermediate bands, 
and -f/sSph includes eight intermediate bands). Calcu- 
lations that restrict this sum are secondarily labeled to 
reflect the restriction. The label "H 14 , no uc" uses -H14, 
but does not include uc bands as intermediate states. 
The label "-H14, no uc/so" uses i/14, but includes nei- 
ther uc nor so bands as intermediate states. The label 
"i?l4, 2BT" uses H14, but only includes two-band terms 
(terms for which the intermediate band is the same as the 
initial or final band). Similar labels are used for -ffssph, 
for example, "TJsSph-PBA, no so" uses the perturbative 
solution to -ffggph and does not include so intermediate 
states. 




100 200 300 400 500 100 200 300 400 500 

Excess photon energy (meV) 

FIG. 2: (color online): Spectra of t)bi (black lines), r)B2 (red 
lines), and rjc (blue lines) for GaAs. Panel (a) shows the 
contributions from each initial valence band; dashed, dotted, 
and dashed-dotted lines include only transitions from the hh, 
Ih, and so bands respectively, while the solid lines include all 
three transitions. The thin solid, light brown line in (a) is the 
total Re (?7?n att ) • Panel (b) separates the total into electron 
(dashed) and hole (dotted) contributions. 



III. CURRENT 

The current injection rate due to the field can be 
written 

T = r$E&E& + r {1) + ^E^E^E^E™, (13) 
where J is the macroscopic current density, and 

jh ) =^!)El*E^E l 2u} + c.c. (14) 

The third rank tensor rf^ k describes one-photon current 
injection (the circular photogalvanic effect) the fifth 
rank tensor T)%~. m describes two-photon current injec- 
tion, and the fourth rank tensor rj&s describes "1+2" 

current injection^ Aversa and Sipe showed that rfA-, is 
related to a doubly divergent part of the third-order non- 
linear susceptibility \^ I n cubic materials with point 
group symmetry Td, Oh or O, a general fourth rank ten- 
sor has four independent components, but due to the 
intrinsic symmetry ijiJ = rfA^ , has only three inde- 
pendent components; there are 21 non-zero components 
of 77(j) in the standard cubic basis: rj^ a = i] bbhb = rj^ c , 

^baab ~abba ~caac ^acca ~cbbc ^beeb 

a(ab)b b(bc)c c(ca)a a(ac)c c(cb)b b(ba)a 

V(j) = -\r> - \r> ~ \r> - \i) 

(the components in parentheses can be exchanged), 
where a, b, and c denote components along the princi- 
pal cubic axesi This can be written 

rjffi = (5 ij S kl + S ik P l ) + ir]B25 u V k + ir] C S m , 

(15) 

where 6 lJ is a Kronecker delta and the only non-isotropic 
part is S 1 ^ , which we define in the principal cubic basis 
as S n 1 = 1 when % = j = k = I and zero otherwise. 



The three independent components are t\b\ = —2irj , 
r]B2 = -irj abba , and ri C = 2irj aabb + irj abba - iif aaa . Thus, 
in a cubic material, 

4) =i VB1 (e: • e 2w ) si* + i VB2 (e u • E^r el 

+ i Vc S ijkl Ei*E^E l 2u +c.c. 

This generalizes the notation we used previously for a 
calculation in the parabolic-band approximation, 3 with 
the connection r\B\ = eDB±/h, and r/B2 = eDB2/h. In 
that, or any other spherical approximation, rjc = 0. 

To calculate f]m, we use the microscopic expression 
first given by Atanasov et al.^ modified to explicitly in- 
clude the sum over spin statesii2i2& An alternate micro- 
scopic expression has been derived in the length gauge, 67 
but it has not yet been used in a calculation. In the in- 
dependent particle approximation that we employ here, 
77(7) is purely imaginary 1 and hence rjsi, i]B2, and rjc are 
real, although they can be complex if excitonic effects are 
included^ - 

The spectra of tjbi, 77/32, and rjc, calculated for GaAs, 
are shown in Fig. EJa) along with the contributions to 
each tensor component from each possible initial valence 
band. For a given photon energy, electrons photoexcited 
from the hh band have higher energies and velocities than 
electrons photoexcited from the Ih band; hence the dom- 
inant component t]bi is larger for hh-c transitions than 
Ih-c transitions. The smallness of r]B2 is due to contri- 
butions from the hh-c transitions having opposite sign to 
the Ih-c transitions, as shown previously in the PBA. 3 

Figure Hfb) separates each tensor component into an 
electron contribution and a hole contribution (denoted 
% and r)h by Atanasov et al). 1 Electrons make a larger 
contribution to r\B\ than holes, due to the lower effective 
mass (and hence higher velocity) of an electron than of 
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a hole (much lower, in the case of a heavy hole) with the 
same crystal momentum. Holes dominate r)B2 at lower 
photon energies, while electrons dominate rjB2 at higher 
energies. Both electrons and holes contribute equally to 
the anisotropic component r\c- 

To help in understanding the importance of the various 
intermediate states, in Fig. [3] we compare the calculated 
current injection tensor elements with various degrees of 
approximation described in Sec. Ill Dl 

The component i}b\ (and hence r]^ aa , since tjbi is 
larger than r]B2 + Vc) is dominated by two-band terms. 
Three-band terms cause an increase, by as much as 34%, 
of t]bi [the difference between the dashed and solid black 
lines in Fig. Ola)]. Although not shown in Fig. |3[ most 
of the increase is due to three-band terms with the so 
band as an intermediate state. Terms with the uc bands 
as intermediate states only cause a small increase to t\b i 
(the difference between the dotted and solid black lines) . 
The warping of the bands is clearly not important for 
TjBi, since the calculation with Hg,s P h closely approxi- 
mates the calculation "i?i4, no uc", which includes the 
same intermediate states. Surprisingly, the "i/gSph-PBA, 
2BT" resulli&^i closely approximates the complete, non- 
perturbative fourteen-band calculation, even at excess 
photon energies for which band nonparabolicity is signif- 
icant. This is due to a fortuitous compensation between 
the neglect of nonparabolicity and the neglect of three- 
band terms. The compensation is not as complete for all 
materials. 

The component r\B2i which determines the current due 
to orthogonal linearly polarized fields, is less forgiving to 
approximations than the component r\B\- We have al- 
ready seen in Fig. [3 that r\B2 is small due to a near can- 
cellation of hh and Ih initial states. Reasonable accuracy 
on r\B2 thus requires higher accuracy on the contribution 
from each initial state. In particular, three-band terms 
must not be neglected. By comparing the dashed-dotted 
and solid lines in Fig.[3fb), it can be seen that, whereas 
the sum of the two-band terms is negative, the sum of 
the three-band terms is positive and of the same mag- 
nitude. It is useful to divide the three-band terms into 
three groups: those with intermediate state from the hh 
or Ih bands, those with intermediate state from the so 
band, and those with intermediate state from one of the 
uc bands. We find that each group contributes roughly 
the same positive amount to r\B2 for excess photon ener- 
gies less than A . The groups are added successively to 
the 2BTs in the dashed, dotted, and solid lines in Fig. 
Efb). Three-band terms with so intermediate states are 
less important at the higher excess photon energies in 
Fig. OJb) . The warping of the bands makes a small but 
non-negligible contribution to r]B2, as seen in the differ- 
ence between the dashed-double-dotted and dotted lines 
of Fig. IHJb). The solid brown line in Fig. EJb) is the 
"#8Sph-PBA, no so" result. 3 At low excess photon ener- 
gies, it greatly underestimates r]B2 due to the neglect of 
so and uc intermediate states, while at excess photon en- 
ergies greater than 100 meV, this is partly compensated 



for by the neglect of nonparabolicity. It appears from 
the difference between "ifgSph-PBA, no so" and "if 14, 
no uc/ so" in Fig.^b) that nonparabolicity becomes im- 
portant at energies above 70 meV. 

The term r\c is purely due to cubic anisotropy by 
definition; in any model that is spherically symmetric 
it is identically zero. There is no cubic anisotropy in 
the "bare" (i.e. without remote band effects) eight-band 
model on the set {r6 C , Tg, v , Tj v }. Cubic anisotropy in 
the fourteen-band model is due to the momentum ma- 
trix elements governed by the parameters Eq and Epi, 
the interband spin-orbit coupling A - , and remote bands 
through (72 —73) and Ck- From Fig.^c), it can be seen 
that three-band terms are important for r\c ■ In fact, with 
only 2BTs included, r/c is positive for GaAs, whereas it 
is negative with all terms included. From Fig.^c) it can 
also be seen that the so band and uc bands are important 
as intermediate states for r\c ■ 

Our calculation of is of the same order of magni- 
tude as the ab initio calculation of Atanasov et al.fi but 
its spectral dependence is different. In particular, tjbi 
agrees more closely with the PBA calculation, as seen in 
Fig-EHa). Atanasov et al. had attributed the difference 
between their ab initio and PBA calculations to the as- 
sumption of k-independent velocity matrix elements in 
the PBA. 1 However, our calculation accounts for the k- 
dependence of velocity matrix elements and agrees closely 
(for rjBi and Re77 aaaa ) to the PBA. The earlier ab initio 
calculation! was, in fact, inaccurate at low photon ener- 
gies due to various computational issues; an improved ab 
initio calculation agrees with the spectral dependence at 
low photon energy given hereS 

Figure 01 shows the spectra of rjBi, r]B2, and r\c cal- 
culated with 7Ji4 for InSb, GaSb, InP, and ZnSe. The 
dashed black line in Fig. El is the PBA resultA^ The 
PBA appears to be a reasonable approximation to rjBi 
for excess energies less than about 0.2E g . In each mate- 
rial, r\B2 <C f]B\-, and in each material except for ZnSe, 
the sign of T)B2 varies as a function of frequency. The 
component r\c-, which arises due to cubic anisotropy, is 
negative for each material. 

The cubic anisotropy of current injection due to co- 
linearly polarized fields can be significant enough that 
it should be measurable. For fields colinearly polarized 
along e, specified by polar angles 9 and <fi relative to the 
cubic axes, 



• e = 21m (-E^oJ Vbi+Vb2 + Vc 



> 



where / (9, 4>) = sin 2 (29) + sin 4 (9) sin 2 (2(f). In general, 
J(j) also has a component perpendicular to e that is pro- 
portional to tjc, but it vanishes for e parallel to (001), 
(110), (111). The field polarization that maximizes the 
current injection depends on the relative sign of rjc and 
Re?7 aaaa = riBi + VB2 + Vc- When they have the op- 
posite sign, current injection is a minimum for e || (001) 
(/ = 0) and a maximum for e || (111) (/ = 4/3); for light 
normally incident on a {001} surface, the largest cur- 
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FIG. 3: (color online): Approximations for GaAs current injection tensor components (a) Wbi, (b) r)B2, and (c) r\c- The 
approximations are described in Sec. Ill Dl 
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FIG. 4: (color online): t]bi (black), r\Bi (red), and r\c (blue) for (a) InSb, (b) GaSb, (c) InP, and (d) ZnSe. The solid lines are 
calculated with the complete fourteen-band model. The dashed line for t]bi is "Hgsph-PBA, 2BT"i2*— • The inset of panel (c) 
shows the area near the origin in more detail. 



rent injection occurs when e || (110) (/ = 1). When they 
have the same sign, current-injection is a maximum for 
e || (001) and a minimum for e || (111). From the GaAs 
results shown in Fig. Ufa), the current injection for the 
three cases e || (001), e || (110), and e || (111) are in the 
ratio 1 to 1.14 to 1.20 at the band edge, 1 to 1.15 to 1.20 
at 200 meV excess photon energy, and 1 to 1.22 to 1.29 
at 500 meV excess photon energy. In contrast, the ab 
initio calculation of Atanasov et al. yields larger ratios, 
for example 1 to 1.32 to 1.43 at 300 meV excess photon 
energy A This disagreement is consistent with the inaccu- 
racy of the ab initio calculation discussed above. Initial 
experiments with GaAs used e || [001] whereas Roos et 
al. exploited the larger signal for e || [HO]^ For each of 
the materials shown in Fig. 01 the minimum current injec- 
tion is for e || (001). It is worth noting that two-photon 
absorption is also a minimum with e || (001) for many 
semiconductors. 69,7 2i2i It seems that both "1+2" current 



injection and two-photon absorption with linearly polar- 
ized fields are larger for e directed along the bonds. 

The cubic anisotropy of "1+2" current injection is pro- 
nounced for cross-linearly polarized fields and opposite- 
circularly polarized fields. For example, for cross-linearly 
polarized fields normally incident on (001) with e w = 
a cos cj) + b sin <fi and §2^ = — a sin <f> + b cos (j>, 



j (J) =lm{ElE^) 



(2r]B2 + Vc sin 2 (2(f>)) e 2a) - -y sin (40) e u 

(18) 

For fields with opposite circular polarizations, the cur- 
rent injection is proportional to r\c an d is hence purely 
anisotropic. 

The component r\c causes a type of current injection 
that has not previously been noted. In all "1+2" exper- 
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iments considered thus far with light normally incident 
on a surface, the direction of current injection lies in the 
plane of the surface. However, with co-linearly polarized 
light fields normally incident on a (111) surface, the cur- 
rent can have a component into (or out of) the surface. 
The current in this case is 



J (I) =2Im(£££ 2 *J 



V2 



j)c cos (38) z 



(19) 



where fj = (j]bi + r\B2 + hvc), z is the [111] direction, 
and 8 is the angle between e and the [211] direction. 
Thus, rjc governs this "surfacing" current. 



IV. POPULATION CONTROL 



The carrier injection rate due to the field can be 
written N = Nm + Nm + ^(2)> where N is the density 
of electron-hole pairs, iVm = ^/LE^E^ 1S one-photon 
absorption, AT (2) = $£ 1 E%E£E*EI, is two-photon ab- 
sorption, and 



%) =Zf I k ) E*jE*jEl + c.c. 



(20) 



i J k 



is "1+2" population control^ The third-rank tensor £ 

has intrinsic symmetry ^jjs — Cm- I n centrosymmet- 
ric materials, such as those with the diamond structure 
(point group Oh), £^ is identically zero; hence, popula- 
tion control requires a noncentrosymmetric material. In 
a material with zinc-blende symmetry (point group T4), 
CI 1 ? has only one independent component; in the stan- 



Ccab tbca each tbac £cba 

Hi) - Hi) ~ HI) ~ HD ~ HD 



Hi) 

dard cubic basis, — — — <, m — sm — ^ m 
are the only non-zero components, where a, b, and c de- 
note components along the principal cubic axes. 

We calculate £(/) with the microscopic expression given 
by Fraser et al., which was derived in the independent- 
particle approximation, and is restricted to Hui < E g < 
2hw& Under those conditions, £m is real and is propor- 
tional to the imaginary part of the susceptibility for sec- 
ond harmonic generation (SHG) 2 ^; specifically, (in mks) 



tabc 

Hi) 



2e 



Imx (2)cba (-2w;w,w). 



(21) 



This connection to SHG, which can be derived from 
considerations of energy transfer and macroscopic 
electrodynamicSf^iiS is important because the imagi- 
nary part of (— 2cj;cj,cj) has sometimes been pre- 
sented en route to a calculation of | .SiZ^SiZSiZLSiZ 9 . 
As well, analytic expressions have been derived for 
the dispersion of SHG by using simple band models, 
with approximations appropriate for 2hw near the band 
gapi£iSL2ii22iSi§i However, these earlier works did not 
connect Im^ 2 -' (— 2uj;uj,lj) with population control, and 



in fact typically stated that it was not independently ob- 
servable. 

Fig. [3] shows the calculation of Imx( 2 ) cfca (— 2u; u>, u) 
for InSb, GaSb, InP, GaAs and ZnSe. Also shown for 
comparison is the PBA expression (|B3(I . derived in Ap- 
pendix [5] Each spectrum can be divided into roughly 
three regions. At very low excess photon energies, visible 
in the log- log plot Fig.^f), the spectrum is roughly inde- 
pendent of uj. This flat part of the spectrum disappears 
if Ck is set to zero; hence, it is due to the k-linear term in 
the c band spin-splitting. Next higher in photon energy, 
up to about 100 meV in GaSb, InP, GaAs, and ZnSe (up 
to about 15 meV in InSb), is a region where the agree- 
ment with the analytic expression l|B3|) is best. In this 
region, the ratio X2/X1, defined in Appendix iBl is 0.37 
for InSb, 0.30 for GaSb, -0.25 for InP, 0.08 for GaAs, 
and 0.07 for ZnSe. At higher photon energies, the dis- 
persion of Imx (2)cfca (-2w; u>,uj) deviates from the PBA 
expression due to band nonparabolicity and warping, k- 
dependence of matrix elements, and transitions from the 
split-off band, which are not included in (|B3(I . 

If we remove the two-band transitions hh-{hh, c}-c, 
lh-{lh, c}-c, and so-{so. c}-c, then the calculation of 
j m ^( 2 ) ( or i s unchanged. This is expected for 

materials of zinc-blende symmetryi2ii2£ Further, many 
years ago Aspnes argued that the so-called "virtual hole 
terms" of the form lh-{so, hh}-c and hh-{so, lh}-c make 
only a small contribution to x^ 2 ) (O)p^ Such terms have 
been neglected in some previous calculations of x^ 
dispersion^ 4 *^ By removing the virtual hole terms, leav- 
ing only {so, Ih, hh}-uc-c transitions, we find £(n is re- 
duced by only 6-10% over the range from the band edge 
to 500 meV above the gap for GaAs. It is thus clear that 
inclusion of the uc bands is necessary for a calculation 
of population control. For some purposes it is also suffi- 
cient, since if remote band effects arc removed from the 
model, leaving the "bare" fourteen-band model ; 44 ' 49 £m 
is decreased by only 7-10% from its full value for GaAs. 

For most materials, the results in Fig. [S] are in 
reasonable agreement with previous calculations of 
I m ^( 2 ) f 73 i 76 i 78 J 79 although most previous calculations had 
poor spectral resolution in this energy range. However, 
for ZnSe, the situation is more complicated. The calcula- 
tion of Huang and Chin is about an order of magnitude 
smaller than ourSfSi and that of Ghahramani et al. is 
about 5 times smaller than oursS Note also that Huang 
and Chin calculated x^ (0) f° r ZnSe to be an order of 
magnitude smaller than experimental results^ Wagner 
et al. have measured the dispersion of |x^l; which is an 
upper bound on Imx^ 2 ^; for ZnSe it is about a factor of 
two smaller than our calculation of Im;^ 2 ) £!i Note that 
Wagner et al. give a different set of band parameters than 
we have used here^ 

The magnitude of £°^ c determines the magnitude of 
population control, but in an experiment one is more in- 
terested in the depth of the phase-dependent modulation 
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FIG. 5: Imx (2) calculated with H14 (solid line) and H14-PBA (dashed line) for (a) InSb, (b) GaSb, (c) InP, (d) GaAs, and (e) 
ZnSe. Panel (f) shows the GaAs calculations on a log-log plot. 



of the carrier absorption, i.e. the control ratio R£ It is 

N(i) Zff ) E»EiEZ+c.c. 
N(i) + N(2) ~ ^K&Kl + ($E*Ei*E*El ' 

This ratio is largest for field amplitudes that equalize 
Nm\ and iV( 2 ^S; in what follows, we assume this con- 
dition has been met. The ratio then depends only on 
£m j £m> C(2) 1 ' anc ^ the polarizations of the two fields. 
For light normally incident on a (111) surface, linearly- 
polarized fields yield R = V2£?if / ^3£ffiff ) aa (1 - a/2), 
while opposite circularly-polarized fields yield 

i? = 2£m/y%)£ ( a 2 7 a (1- a/6-5), (22) 

where a = (£ffi a - ^ b - 2^ b )/^ ) aa and 6 = 
^aaaa + £aabb _ 2^ b ) / (2£™ ) aa ) are two-photon absorp- 
tion anisotropy and circular dichroism parameters'^ 
Stevens et al. found that for light normally incident on a 
(111) surface of GaAs, opposite circularly polarized fields 
yield the largest ratioASl For light normally incident on 
a (110) surface, fields linearly polarized along [111] yield 

R = 2£ m /^3£m£f 2 a ) aa (l-2<7/3). (23) 

The polarization configuration that yields a global max- 
imum for the control ratio depends on the material and 
photon energy; we have found that (1221) is the maximum 
except for very close to the band edge, where 12.'iH is the 
maximum. 

To calculate the population control ratio, it is desir- 
able to use values of £m, £m, and £(2) calculated within 
the same set of approximations. We use microscopic 
expressions for £(1) and £(2) in the independent-particle 



approximation^ 1 and calculate them within the fourteen- 
band model. Note that our calculation of two-photon 
absorption (£(2)) is similar to that of Hutchings and 
Wherrett^ but that our model includes remote band 
effects. 

Fig. shows the calculated spectra of the population 
control ratio (|22|l for various semiconductors. For each 
material, the ratio is close to unity at the band edge, 
then drops steeply, but flattens out to some non-zero ra- 
tio as photon energy is increased. In general, the smaller 
the band gap (or conduction band effective mass) of the 
material, the narrower the range over which the ratio 
drops, and the lower the ratio at higher excess photon 
energy. Worth noting is the particularly large ratio for 
ZnSe. Also plotted in Fig. [S] is the ratio appropriate 
for linearly-polarized fields normally incident on a (111) 
surface of GaAs, which was the configuration in the ex- 
periment of Fraser et ali£ For all materials, the ratio 12.')!) 
reaches exactly unity at the band edge, in agreement with 
the PBA calculation (jB9|) in Appendix |B| 

The only previous theoretical calculation of the 
population-control ratio, which was for GaAs, missed 
finding the large ratio near the band edge because it was 
based on ab initio calculations of £m, £(2) and that 
had poor spectral resolution near the band edge^ Over 
the rest of the spectrum shown in Fig. EI it is about a fac- 
tor of two smaller than our calculation. This is consistent 
with the previous calculation being based on a calcula- 
tion of the two-photon absorption coefficient £(2) that is 
too large by comparison with other calculations. ^'71 

The population-control ratio has been measured only 
in GaAsi 2 i 4 i 9 i 27 i 30 The measured ratios on (lll)-GaAs, at 
excess photon energies of 180 meV^S and 312 meV^Si 
were 4 to 5 times smaller than our calculation. Some 
of the difference can be attributed to phase mismatch 
and large sample thickness £2i2LSI An experiment on a 
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For a cubic material, \i^ m nas 54 non-zero elements in 
the principal cubic basis, and can be written 

ijklm A*JV1 / jmlcik , _j'mfc cil\ , ,, Ajm\kl 
[ifa = — — (S J 6 + S J ) + a 
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Excess photon energy (meV) 

FIG. 6: (color online) Calculated population control ratios 
appropriate for opposite circularly polarized fields normally 
incident on a (111) surface of InSb, GaSb, InP, GaAs, and 
ZnSe. The blue, dotted line is the ratio for linearly polarized 
fields normally incident on a (111) surface of GaAs. 



MW2 



( £ iml 5 jk + £ imk 8 jl^ + ^^ikln £ njm 



2 (26) 

where the non-isotropic tensor 8 l ^ kl has nonzero compo- 
nents 5 aaaa = 8 bbbb = 8 CCCC = 1, where a, 6, and c denote 
components along the principal cubic axes. The six inde- 
pendent components are /ijvi = 2/^ aba , ^N2 = 2/z^ aba , 

I-IN3 = Mf/) CCC ! MCI = H(I) aC - UNI ~ HN3, HC2 = 
pbaaac _ ^ + and ^ = 2^f ) acb - fl N1 - fl N2 . 

Thus in a cubic material, 



(llO)-grown multiple quantum well was complicated by 
an additional cascaded second harmonic effect £ 



V. SPIN CURRENT 



K\] =vniE*J (E 2w x Elf + m2 (E 2w x E*J 



+ c. 

(27) 



_i_ I,, xikln njrn . xjkln nim\ Tp*k rp*l rpm ■ 



Spin-current density can be quantified by a second- 
rank pseudotensor lf y defined as the average value of 
the product v l S : > , where v is the velocity operator and 
S is the spin operator^ Note that some authors alter- 
nately choose the first index to represent spin and the 
second index to represent velocity. 87 Also, due to the 
spin-orbit part of the velocity operator — the so-called 
"anomalous" velocity [the second term in J2J] — v and 
S do not commute, and thus v l S J is not Hermitian. In- 
stead, one should take (v l S 3 + Siv l )/2 as the operator 
for spin-current. But since we neglect the anomalous ve- 
locity (see Appendix |A"|) . this is not necessary. 

The spin-current injection rate due to the field (IJ can 
be written 



ijkl 



i .Ajklmn jpk* jpl* jprn ipn (oa\ 



where the pseudotensor describes one-photon spin- 
current injection^ the pseudotensor fi%, 
two-photon spin-current injection, and 



ijklran 



describes 



A (7) 



ijklm jr*kj?*l j?m 
' (Z) W '■d 2lj 



c.c. 



(25) 



Note that the injection of (v ■ S) is zero in a cubic 
material, i.e., K lJ is traceless. In an isotropic model, 
such as the one we used previously^ /ici = fJ-C2 — 
MC3 = 0. The connection to our previous notation is 
Mjvi = D(A\-Ai), fi N2 = D(A 2 + A i ), and ^ 3 = 
D(A 3 + A 4 )£ 

The spin-current injection can be divided into a con- 

and a contribution from 



tribution from electrons KVj. 
holes K%. h y that is, K%s = 

ijklm ijklm . ijklm 



(I) 



fl 



(I;h) 



)■ 



(similarly, 
Expressions in the PBA 



for both the electron and hole spin current are given 
elsewhere 3 -; here we focus on the electron spin current, 
since hole spin relaxation is typically very fastiS2i2S 

A microscopic expression for the spin-current injection 
was derived previously in the Fermi's golden rule (FGR) 
limit of perturbation theory and applied to a model in 
which all bands are doubly degenerated However, it is 
unsuitable for a calculation with if 14, which accounts for 
the small splitting of the spin degeneracy that occurs 
in materials of zinc-blende symmetry £££±£1 If the spin- 
split bands were well separated, then the microscopic ex- 
pression for -^(/-e) would be 



is "1+2" spin-current injections 



1 , ijklm 1 1 ijlkm 

dotensor /j,/^ has intrinsic symmetry 

t • • t - 1 ijklm 

In an isotropic material, 



The fifth-rank pseu- 

ijklm 

(i) ~ • 
T j has three independent 
components, while in a cubic material (with T^, O, or Oh 
symmetry) /x!^ im has six independent components. The 
four parameters Ai, i = 1-4, that we used previously to 
describe spin-current injection in an isotropic model 3 can 
be reduced to three independent components with iden- 
tities such as e l 3 m 5 kl - s ijk 5 lm + £ j km S il - £ lkm S jl = 0^ 



(I;e 



2tt 

= L3 



<ck| v i S i |ck) 



5 (2oj - u cv (k)) , 



where L 3 is a normalization volume; the one-photon am- 
plitude Q^l k is 



^ = l 2^ E2 - V - (k) - 



(28) 
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where the charge on an electron is e (e < 0), and the 



(2) 

two-photon amplitude Q, c v k is 



(2) 



( 



e 



(E h ,-v C|n (k))(E^-v n ,„(k)) 
w nv (k) - u 



(29) 



However, for the photon energies and materials studied 
here, the spin-splitting is small; it is comparable to the 
broadening that one would calculate from the scattering 
time of the states, and also to the laser bandwidth for 
typical ultrafast experiments. Thus, the spin-split bands 
should be treated as quasidegenerate in FGR, with the 
result 



k ij =— 



i 



X - [6 (2uj - uj cv (k)) + 5 (2uj - uj dv (k))] + c.c., 

where the prime on the summation indicates a restric- 
tion to pairs (c, d) for which either d = c, or c and d 
are a quasidegenerate pair. The optical excitation of the 
coherence between spin-split bands can be justified using 
the semiconductor optical Bloch equation approach, as 
was done for the one-photon spin properties^ Note that 
this issue does not arise for "1+2" current injection or 
"1+2" population control, since (ck|v|c'k) and (ck|c'k) 
vanish between spin-split bands. 

Using the time-reversal properties of the Bloch func- 
tions, we find that /U(j ;e ) is real, and can be written as 



c,c' v,k n 

(ck\v l S^ |c'k) 



x Re 



[i\,j-klm { -\jhlra \* 

c,c\v V c'.c.v) 



(30) 



where 



klm 
c.c' .V 



M, 



-v? v (k) [v k c *(k)vt(k)+v l c * n (k)vt(k)] 



(31) 

That /i(/ ;e ) in (|30[1 is purely real is a consequence of the 
independent-particle approximation. 40 



A. Calculation results 

The spectra of the independent components of /X(/ ;e ), 
calculated for GaAs, are shown in Fig. [7] and Fig.[S] Fig- 
ure \7\ also shows contributions from each possible ini- 
tial valence band. Figure [S] shows the spin-current injec- 
tion calculated with various degrees of approximation de- 
scribed in Sec. Ill Dl The only other calculation of "1+2" 
spin-current injection for bulk GaAs is our earlier calcu- 
lation, which used a spherical, parabolic-band approxi- 
mation to the eight-band model and did not include the 
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FIG. 7: (color online): Calculated spectra of GaAs spin- 
current injection components and their contributions from 
each initial valence band; dashed, dotted, and dashed-dotted 
lines include only transitions from the hh, Ih, and so bands 
respectively, while the solid lines include all three transitions. 
Panel (a) shows /iati (black lines), /ijv2 (red lines), and fiN3 
(blue lines). Panel (b) shows fici (black lines), fic2 (red 
lines), and fic3 (blue lines). 



so band as an intermediate state- ; it is shown in Fig. [S] 
for tx N1 , fi N2 , and fi N3 . 

The term fi^i has the largest magnitude of the six in- 
dependent parameters of H(i- e ) ■ Since it is negative for 
hh and Ih transitions but positive for so transitions, it 
peaks in magnitude at 2hu just above E g + An. (the en- 
ergy at which so transitions become allowed) . Two band 
terms make the largest contribution to /Ltjvi, followed by 
three-band terms with hh or Ih intermediate states. The 
so and uc intermediate states make a very small contri- 
bution to fj.Ni for excess energies less than 200 meV. The 
warping of the bands is not important for (J,ni, since the 
calculation with -ffsSph closely approximates the "i?i4, 
no uc" calculation, which includes the same intermediate 
states. The "TJsSph-PBA, no so" calculation, which we 
derived previously^ is a good approximation to /ijvi at 
excess energies below 250 meV; nonparabolicity becomes 
important at higher energies. The hh contribution has 
a larger magnitude than the Ih contribution in part be- 
cause three-band terms increase the magnitude of the hh 
contribution, but decrease that of the Ih contribution, as 
expected from the PBA expression (IBla|) . 

The term fi^2 is negative for hh transitions, positive 
for Ih transitions, and negligible for so transitions. The 
calculation "i?i4, 2BT" is a good approximation to the 
calculation H14. However, the three-band terms are not 
small; rather, they nearly cancel. In particular the tran- 
sition hh-lh-c makes a large positive contribution to HN2, 
while the transition hh-so-c makes a large negative con- 
tribution. Since our earlier PBA calculation included the 
former but not the latter^ it is a poor approximation 
to /J-N2- But by including only 2BTs, it is a fair ap- 
proximation for excess energies less than 200 meV. This 
agreement is fortuitous, since the calculation -f/ssph un- 
derestimates the magnitude of fiN2, and the PBA leads 
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FIG. 8: (color online): Approximations for GaAs spin current components (a) uni, (b) fiN2, (c) /ijv3, (d) fioi, (e) (f) 
HC3- The approximations are described in Sec. Ill Dl 



to an overestimation of the magnitude of /ijV2- 

The term fi^3 is negligible when only 2BTs are in- 
cluded, in agreement with the PBA£ The hh-lh-c transi- 
tions are positive, while the Ih-hh-c transitions are nega- 
tive; the former is larger, and thus is positive when 
so intermediate states are neglected. Both Ih-so-c and 
hh-so-c are negative and substantial enough to make 
the total /iN3 negative. Consequently, our earlier PBA 
result,^ which neglects so intermediate states, is a poor 
approximation to fiN3- Upper conduction bands make a 
fairly small contribution to /iN3, and warping does not 
seem to be important for /xjv3 since the calculation with 
#8Sph is a good approximation. 

As expected, the terms nci, and nc3 are zero 

when calculated with -f/gsph- 

The term fici is negligible when only 2BTs are in- 
cluded. Transitions with intermediate states in the set 
{hh,lh, so} comprise roughly two-thirds of [ici- The 
anisotropy of these transitions is not simply due to the 
warping of the hh and Ih bands, which we have deter- 
mined by a calculation (not shown) using Hg without 



the remote band contribution to the velocity. Rather, 
it comes from wave function mixing of the Tg, c and T7 C 
states into the valence and c band states. The cubic 
anisotropy of two-photon absorption has been attributed 
to such wave function mixing42<S2 The other third of the 
full nci is due to transitions with the uc intermediate 
state, which would be forbidden close to the T point if 
the material were isotropic. We also note that each three- 
band term makes a positive contribution to /ici- 

The term \xci is nearly negligible when only 2BTs are 
included. Transitions from the hh and Ih bands have op- 
posite sign, and those from the so band are negligible. 
About half of fic2 is due to the transitions hh-lh-c and 
Ih-hh-c, and the other half is due to transitions with the 
uc intermediate states. Transitions with so intermediate 
states are negligible. As with /xci, the anisotropy of the 
hh-lh-c and Ih-hh-c transitions is due to the wave func- 
tion mixing of the T% c and Ti c states into the hh, Ih, and 
c band states. 

The term fic3 is positive for hh transitions, negative 
for Ih transitions, and negligible for so transitions. The 
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transitions hh-so-c and Ih-so-c account for most of the 
value of /ic3, but 2BTs are not negligible. Transitions 
with uc intermediate states reduce the value of nc3 by 
as much as 10%. Most of /zc3, especially at energies less 
than 200 meV, is due to the warping of the hh and Ih 
bands. Consistent with this, we find that remote band 
effects are somewhat important for yUc3! when remote 
band effects are removed, the calculation of [ic3 is about 
25% larger than the full calculation. Note that [ic3 is 
far more sensitive to remote band effects than any other 
optical property calculated in this article. 

In Fig.[5]we plot the spectra of the independent compo- 
nents of the spin current density pseudotensor for InSb, 
GaSb, InP, and ZnSe. The spin current tensor is largest 
for InSb in agreement with the PBA expressions in Ap- 
pendix^] We also note that /j,n3 is positive for InSb and 
GaSb at low excess photon energy, whereas it is negative 
for InP, GaAs, and ZnSe. 



B. Configurations 

Co-circularly polarized fields generate a spin-polarized 
current, which can be characterized by its degree of spin 

polarization / = (2e/K) Kl 3 cfW /\J e \, where n is a unit 
vector normal to the polarization plane of the fields, and 
q is a unit vector in the direction of J e £ Essentially, 
/ = (vS) I (v). Since this measure aims to character- 
ize the photoexcited distribution of electrons, we neglect 
holes from both K and J in this calculation^ For fields 
normally incident on a (001) surface (i.e. e u = e 2u = 



(x ± iy) /V2), the spin current is 



Uni H ^~ 



rh l ± z J + I \iN2 



2 ) 



where m± = sin (2</> w - 4> 2uJ ) x ± cos (2(f> w - 4>2u)y, the 
current is J(j) = y/2E 2LlJ E 2 (r]Bi + rjc /2) m±, and the 
degree of spin polarization is 



/ = 



2e /j An + mct/2 
Vbi + Vc/2 ' 



(32) 



For fields normally incident on a (111) surface, J(/j 
\f2E 2u El (ijbi + ?7c/3) m±, and 



/ 



2e /XjV! + nci/3 + ^C3/3 
ft Vbi+Vc/3 



(33) 



The degree of spin polarization is plotted for GaAs in 
Fig. llOf a). The cubic anisotropy is small, but clearly 
seen, especially at low excess photon energies. The other 
materials have very similar degrees of spin polarization. 

A pure spin current, without an electrical current, 
can be generated with cross-linearly polarized fields^ We 
consider fields polarized in the (001) plane, with the oj 
field polarized at an angle 8 to the x axis (i.e. [100]) 
and the 2lo field polarized at an angle 8 to the y axis 
(e w = x cos 8 + y sin 8 and e 2w = — x sin 8 + y cos 8) . The 
spin current is 



K lJ - 



- ^\E 2(Aj \\E^\ 2 cos (2^ - (j) 2uJ ) [(4^vi + 4/^3 + 3/xci + Mci cos (46>))<,P 

- sin (48) (fxcie 2u z j + ^ C 2Z l e J 2l ^j + (4^i N2 ~ A/i N3 + 3^ C 2 + VC2 cos (48)) z l e>^ 



r 



This pure spin current is typically measured by the re- 
sulting displacement of up and down spins££ The fi- 
nite displacement results from transport and scattering 
of the electrons. Using the Boltzmann transport equa- 
tion in the relaxation time approximation with space- 
charge effects justifiably neglected^ one finds cf (z) = 

(4T/h)k^^l (n {1) + Here, d (z) is the dis- 

placement of spins measured with respect to the quan- 
tization direction z, and r is the momentum relaxation 
time. We assume the field intensities have been chosen 
to balance one- and two-photon absorption, a condition 
that is ^-dependent due to the cubic anisotropy of two- 



photon absorption. Thus, 

r (4/xjvi + 4^jv 3 + 3/zci + Mci cos (48)) 



d (z) • 



and 



d(z) • e 2u 



h ^^(1- (a/2) sin 2 (28)) 
t nci sin (48) 



sft^T (1- (a/2) sin 2 (28)) 



(34) 
, (35) 



where a is the two-photon absorption cubic-anisotropy 
factor given explicitly in the next section42iS£ At = 
and 8 — tt/4, d is parallel to e u . The spin separation 
distance is plotted in Fig. HOf b) , where we have assumed 
a momentum relaxation time of 100 fs for each material. 

This calculation of the spin separation distance is a 
significant improvement over our initial calculations^^ 
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FIG. 9: Spectra of spin current components for InSb, GaSb, InP, and ZnSe: /j,ni (solid black line), fiN2 (solid red line), fiN3 
(solid blue line), fici (dashed black line), fiQ2 (dashed red line), and fic3 (dashed blue line). 




Excess photon energy (meV) 

FIG. 10: (a) Degree of polarization of spin-polarized current 
due to co-circularly polarized fields, (b) Displacement of spins 
in pure spin current due to cross- linearly polarized fields. 

which used the eight-band PBA and neglected three- 
band terms from the two-photon amplitude ("-ffssph- 
PBA, 2BT"). Stevens et al. measured a spin separation 
distance of 20 nm in a GaAs multiple quantum well at an 
excess photon energy of 200 meV, and estimated a mo- 
mentum relaxation time of r = 45 fs. 7 For r = 45 fs, we 
calculate a spin separation distance of 20.0 nm for bulk 
GaAs at 200 meV. Hiibner et al. measured a spin sep- 
aration distance of 24 nm (the photoluminescence spot 
separation is half this distance) in cubic ZnSe at an excess 
photon energy of 280 meV, and estimated a momentum 
relaxation time of r = 100 fsp2> The calculation in Fig. 
ITUTb) yields d = 23.6 nm for ZnSe at 280 meV. In both 
cases, we now find very good agreement with the experi- 
ment, whereas the previous model resulted in larger spin 
separation distances. Of course, this agreement is contin- 
gent on the accuracy of the momentum relaxation time 
estimates. 

Note that both the degree of spin polarization for co- 
circularly polarized fields and the spin-separation dis- 
tance, plotted in Fig. have a kink at excess photon 
energy A and decrease at higher excess photon energies. 
A similar kink and decrease, due to the onset of transi- 



tions from the split-off band, occurs for both one-photon 
spin injection** and two-photon spin injection^ 

VI. SPIN CONTROL 

The spin injection rate due to the field can be 
written S = Sm + Sm + S(2)j where S is the macro- 
scopic spin density, Slg = C~L E^ is one-photon spin 
injection, 94 S\ 2) = C$ lm E% E** E l u E™ is two-photon 
spin injection, 52 and 

SIj^C^K'K'eL+c.c. (36) 

is "1+2" spin control. 4 In previous sections and in some 
of the expressions below, we use S to denote the single- 
particle spin operator. It should be obvious by context 
when S refers to the macroscopic spin density and when 
it refers to that spin operator. 

The fourth-rank pseudotensor £^ has intrinsic sym- 
metry on the indices j <-» k. Such a pseudotensor is zero 
in the presence of inversion symmetry; hence, "1+2" spin 
control requires materials of lower symmetry. For zinc- 
blende symmetry (point group T c i) , a general fourth-rank 
pseudotensor has three independent parameters and 18 
non-zero elements in the standard cubic basis; forcing the 
j <-> k symmetry leaves two independent parameters 

■/- /-abba /-caac /~bccb /-acca /-ebbe /-baab 

l( ,IA = 4(7) - 4(7) - 4(7) - -C,( Z ) - -(,(/) - -(,(/) 

(37a) 

■/• /-aabb /-ccaa /-bbec /-aacc /-cebb /-bbaa 

i(,ib = (,(/) - C,( 7) - - - - 

Cabab /-caca /-bebe /-acac /-cbcb /-baba 

7) - C(J) - 4(7) - -4(7) - -4(7) - -4(7) 

(37b) 

The spin injection has a contribution from electrons 
S(7 ;e ), and a contribution from holes S(7 ; /j); that is, 
S(7) = S(7 ;e ) + S(7 ;/l ), and £(7) = C(7 ;e ) + C(i-,h)- 
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FIG. 11: (color online) Spin control pseudotensor components 
£ja (black lines), (^ib (red lines), and (£ia + %Cib) (blue line) 
with breakdown into initial states. Dotted lines include tran- 
sitions from the Ih band, dashed lines include transitions from 
the hh band, dashed-dotted lines include transitions from the 
so band, and solid lines include all transitions. 

We treat the spin-split bands as quasidegenerate when 
taking the FGR limit of perturbation theory, as discussed 
for the spin current in Section^ deriving the microscopic 
expression 

x-[iS (2w - uj cv (k)) + 5 (2cj - uj c , v (k))] + c.c, 

where the prime on the summation indicates a restric- 
tion to pairs (c, c') for which either d — c, or c and d are 
a quasidegenerate pair. Using the time-reversal proper- 
ties of the Bloch functions, we find that C(J";e) is purely 
imaginary and can be written 

c,c' u,k n 

xR e ((|l^)[ M « ( <; )•]! 

(38) 

where M;? c , ^ is given in Eq. l(3*T|) . 

The spectra of Qa and C/s for GaAs are shown in Figs. 
IlllandlT^l Figure ITTI also shows the contributions from 
each possible initial valence band. Figure IT21 shows the 
spin control calculated with various degrees of approxi- 
mation described in Sec. Ill Dl 

The term £ia decreases from zero at the band edge to a 
maximum negative value at 40 meV, mostly due to tran- 
sitions from the hh band, and is positive at higher excess 
photon energies, mostly due to transitions from the Ih 
band. The low energy behavior is in agreement with the 
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FIG. 12: Spin control pseudotensor for GaAs. The calcula- 
tions are i/14 (solid black lines), u Hi4, no uc" (dotted lines), 
"H\4, no uc/so" (dashed-dotted lines), "-H14, 2BT" (dashed 
lines), and ".H14-PBA" (solid grey lines). 



PBA result l|B2ajl . in which the ratio of hh : Ih tran- 
sitions is (m Cj ^^/m Ci i/j) 3 ^ 2 . Transitions with so and uc 
intermediate states dominate the decrease in £7,4 at low 
excess photon energies, as seen in Fig. I12f a); they are 
the only non-zero transitions in the PBA result l|B2a(l . 
The contribution from uc intermediate states is negative 
and approximately constant over most of the spectrum, 
whereas the contribution from so intermediate states 
changes from negative to positive as transitions from the 
so band become allowed (2hui > E 9 + Aq). The contribu- 
tion from 2BTs, which is zero in the PBA, is positive over 
the whole spectrum. The breakdown of the PBA is due 
to the increase in magnitude of the 2BTs. In fact, the 
sum of the PBA and the 2BTs is a good approximation 
to the full calculation. We also note that a calculation 
with Hg for Qa yields a nearly negligible result; thus, 
the contribution from intermediate states within the set 
{so, Ih, hh, c} (including 2BTs) is due to the mixing of 
the T-j c and Tg, c wavefunctions with these states. 

The term Q b is larger in magnitude than the term Qa 
over most of the calculated spectrum. It falls to a max- 
imum negative value at 95 meV, sharply increases when 
transitions from the so band become allowed, and is pos- 
itive at higher excess photon energy. At lower photon 
energies, transitions from the hh band and transitions 
from the Ih band both make negative contributions to 
Cib', in the PBA result (|B2 b|) the ratio of hh : Ih tran- 
sitions is (m c ,hh/nT-c.ih) 3 ^ 2 ■ Fig. I12f b) reveals that £xb 
is essentially due to contributions from uc intermediate 
states, and 2BTs. Over the whole spectrum, the former 
are negative while the latter are positive. The smallness 
of the contribution from so intermediate states is also 
seen in the PBA result (|B2b|l . since Z + > Z'_ in that 
expression. We also note that a calculation with Hs for 
Cib yields a nearly negligible; thus, the contribution from 
intermediate states within the set {so, Ih, hh, c} (includ- 
ing 2BTs) is due to the mixing of the Tj c and Tg c wave 
functions with these states. 

We have also calculated the spin-control pseudotensor 
for the semiconductors InSb, GaSb, InP, and ZnSe. The 
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FIG. 13: (color online) Spin control calculated for InSb, GaSb, InP, and ZnSe. Black lines are (,ia and red lines are Qib- Solid 
lines are the full calculation using -H14. Dotted lines are the parabolic-band approximation calculated with lB2at and JB2bl . 



results are shown in Fig. ^| along with the parabolic- 
band approximations (|B2a(l and (|B2b(l . 

The magnitude of spin control is determined by Oi)> 
but in an experiment one is more interested in the depth 
of the phase-dependent modulation of the spin polariza- 
tion signal. One possible definition for the signal is the 
ratio of spin injection measured with both u> and 2uj fields 
to the sum of the spin injections measured with circularly 
polarized fields of each frequency^ The amplitude of its 
modulation is 



(39) 



(2) 



where the argument (<j + ) indicates injection with a a + 
polarized field. This ratio, which is largest for field am- 
plitudes that equalize S?^ (c + ) and S* 2 -\ (f + )j wa s mea- 
sured by Stevens et al. with excess photon energies of 150 
meV and 280 mcVi£ 

The ratio (|39[1 has an undesirable feature: it can ex- 
ceed unity. Close to the band edge in many semiconduc- 
tors (at 2 meV in GaAs), there is a photon energy for 
which (er + ) = 0. 52 At that photon energy, it is im- 
possible to choose field amplitudes to balance one- and 
two-photon spin injection with circular polarized fields 
[Sfa (c + ) = S? 2 s (<r + )], and thus the maximum ratio has 

a singularity. Even if the condition S'(i) (er + ) = S( 2 ) (°" + ) 
is relaxed, the ratio (|39[1 can exceed unity. This is be- 
cause Svj, (cr + ) and SLs (er + ) have opposite sign close to 

the band gap^ and thus it is possible, by appropriate 
choice of field amplitudes, to make the denominator of 
the ratio arbitrarily small. 

An alternate ratio to characterize the spin control, 
which has an upper bound of unity, is 



Rs 



qz 



kN (1) +N (2) 



(40) 



It is the amplitude of phase-dependent oscillation of the 
degree of spin polarization, and it is most useful when 
there is little or no population control. We assume the 
fields are chosen to balance one- and two-photon absorp- 
tion. For most photon energies and materials this is 
nearly the same as balancing one- and two-photon spin 
injection. 

For normal incidence on a (111) sample, opposite cir- 
cularly polarized fields yield 



Rs 



2 

Ji 



\(ia + 2(ib\ 



For normal incidence on a (110) sample, opposite circu- 
larly polarized fields 



Rs 



2 

^4 



3|C/a + 2C/b| 



and orthogonal linearly polarized fields (xy-polarized) 
yield 



Rs (a) = r- 



|(C/a + 2C/s) (r + 3sin 2 a) 



cos a 



2 V^i a )^( a 2T a ( X ~ 2" CT ( sin2 ") (! + 3 cos 2 a)) 

where r = —2(ia/((ia + 2C/b), 4 and a is the angle be- 
tween the polarization of the u> field (E w ) and the [001] 
axis, which lies in the (110) plane. The angle that max- 
imizes Rs depends on photon energy through r and a. 
We determine it numerically. 

The ratio Rs for GaAs is plotted in Fig. Hlf a). For 
(lll)-incidcncc, opposite circularly polarized fields yield 
the highest ratio over the studied range of photon en- 
ergies. For (llO)-incidence, opposite circularly polarized 
fields yield the highest ratio, except for between 190 meV 
and 415 meV when rry-polarized fields the highest ratio. 
For xy-polarized fields, the angle that yields the largest 
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FIG. 14: (color online) Spin-control ratio normalized by car- 
rier population [Eq. 140H 1. In (a), for GaAs, black lines are 
(lll)-incident, opposite circularly polarized fields; green lines 
are (llO)-incident, opposite circularly polarized fields; and red 
lines are (llO)-incident, orthogonal linearly polarized fields. 
In (b), for InSb, GaSb, InP, GaAs, and ZnSe, solid lines are 
(lll)-incident, opposite circularly polarized fields. The dot- 
ted line in (b) is (llO)-incident, orthogonal linearly polarized 
fields for InP. 



dieted by the fourteen-band calculation close to the band 
edge, where allowed-allowed two-photon transitions dom- 
inate allowed-forbidden two-photon transitions. Neither 
"1+2" population control, nor "1+2" spin control have 
yet been experimentally studied in that spectral range. 



APPENDIX A: NEGLECT OF THE 
ANOMALOUS VELOCITY AND k-DEPENDENT 
SPIN-ORBIT COUPLING 

The anomalous velocity, i.e. = (v — p/m) = 
H(cr x VV) / (4m 2 c 2 ) , which leads to k-dependent spin- 
orbit coupling in from the term Kk ■ ~va , is often ne- 
glected in k • p models ^USJ^/S^i Some authors have 
treated matrix elements of W as additional independent 
parametersj^ 2 *2ii2L2& For example, Bahder, who gives the 
matrix for fck ■ within the eight-band model, defines 
the model parameterize 



ratio decreases from 0.99 rad to 0.53 rad from the band 
edge to 320 meV, and is zero for higher excess energies. 

The ratio Rs for the five semiconductors InSb, GaSb, 
InP, GaAs, and ZnSe are plotted in Fig. UlTbV At low 
photon energy, opposite circularly polarized fields nor- 
mally incident on (111) yield the largest ratio for InSb, 
GaSb, GaAs, and ZnSe, whereas orthogonal linearly po- 
larized fields normally incident on (110) yield the largest 
ratio for InP. 



VII. SUMMARY 

We have studied the four "1+2" coherent control 
effects — current injection, spin-current injection, popula- 
tion control, and spin control — in bulk semiconductors 
having zinc-blende symmetry. We used an empirical, 
fourteen-band k • p Hamiltonian and examined the rela- 
tive importance to each effect of the possible initial and 
intermediate states. We have also studied the crystal 
orientation and polarization dependencies of each effect. 
Cubic anisotropy is small in some cases, but large in oth- 
ers. 

We have compared the numerical calculation with an- 
alytic expressions, derived in the parabolic-band approx- 
imation, to show the value and limitations of the latter. 
The PBA expressions, where they are accurate, are useful 
to show how the effects scale in different materials. 

The comparison between the two approaches estab- 
lishes that, at low excess photon energies, "1+2" cur- 
rent injection and "1+2" spin-current injection are 
due to interference of allowed one-photon transitions 
and allowed-forbidden two-photon transitions, whereas 
"1+2" population control and "1+2" spin control are 
due to interference of allowed one-photon transitions 
and allowed-allowed two photon transitions. It also ex- 
plains the large population- and spin-control ratios pre- 



C = 



1 



4ra 2 c ; 



(S\V X V\X). 



Ostromek used the value Co = 0.16 eV A to fit the eight- 
band model to experimental results^ We here relate ma- 
trix elements of W (and hence matrix elements of va) 
to other parameters of the model, thereby demonstrating 
that they can be neglected for the effects we consider. 

Bir and Pikus showed that the identity [Hq, p] = itiVV 
leads to {X\ V y V \Z) = 0. 100 An application of that iden- 
tity to the remaining nonzero matrix elements yields 



(S\V X V\X) = ^(E S -E X ), 

m P' 

(S\W x V\x)=--^(E x ~E s ). 



iX\ V„Vh) : : (Z|V„I» '^-{1 - E s ! 



(Ala) 
(Alb) 

(Ale) 



and similar results for cyclic permutations and Hermitian 
conjugates of these. The energies Es, Ex, and E x are the 
eigenvalues of \S), \X), and \x) with respect to the Hamil- 
tonian Hq. Their values are fixed by the requirement that 
the eigenvalues of i?k=o yield the parameters E g , E' , Aq, 
and A' ml Neglecting the small contribution from A - , 
E s - E x = E g + A /3, E x - E s = E' - E g + 2A /3, 
and E X -E X =E' + 2A /3 + A /3. 

Thus, (|Ala(l gives matrix elements of 'VV in terms of 
other model parameters. In particular, with parameters 
from Tabled for GaAs, we find C = 5 x 10~ 6 eVA. 

From the point of view of the theory of 
invariants ^iS&iSLiSl k-dependent spin-orbit coupling 
amounts to using different values of Pq for and T7 
valence bands (and similar changes for Pg coupling and 
Q coupling) In terms of C , P -»• P 7 = P + 2v / 3C 
for couplings with T7 bands and Pq —> Pg = Pq — v3Cq 
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for couplings with Tg bands. From HAlaJl . 

Pr-Ps = V3Co = 3{E S -E X ) ^ 3E g - A 
Pq Pq Amc 2 Amc 2 

This is very small, since the numerator is on the order of 
eV, whereas mc 2 = 5.11 x 10 5 eV. And since this relative 
change in the matrix element depends on the ratio of 
Co to Po, even the overly large coupling value of Co = 
0.16 eVA has only a small effect on optical properties 

For comparison, consider interband spin-orbit coupling 
parameterized by A - . In the eight-band model, inter- 
band spin-orbit coupling is a remote band effect (since 
it is a coupling with the uc bands), which effectively 
causes^ -> P7 = Po + (2A~P^) / [3 {E' + A )] and 
Po -> P 8 = Po ~ (A-Po) / [3 (E' + A )]. Thus, 

P 7 - P 8 _ A- P> 
Po ~ E' Po' 

This effect, which is included in our calculation, is small 
(it is 4 x 10 -3 in GaAs), but it is orders of magnitude 
larger than the relative change due to k-dependent spin- 
orbit coupling. 

The above suggests that k-dependent spin-orbit cou- 
pling can be neglected for the processes we consider in 
bulk, cubic materials. 

APPENDIX B: PARABOLIC BAND 
APPROXIMATIONS 

In this appendix, we discuss parabolic-band approxi- 
mation (PBA) expressions, which are perturbative in the 



Bloch wave vector k, for "1+2" coherent control effects. 



1. Current 



There have been several different calculations of 77 
in the PBAiiS&ffi Using a two-band model (one con- 
duction and one valence band), Atanasov et al. ob- 
tained r/Bi <x (2huj — E g ) 3 ^ 2 and r/B2 — OA Using a 
three-band model, but only accounting for two-band 
terms, Shiek-Bahae studied the approximate scaling of 
"1+2" current injection spectra with the band gap E g 
and concluded that tjbi and r)B2 are proportional to 
E- 2 {2x - 1) 3/2 (2xy 4 , where x = hu/Eg 38 Our earlier 
PBA calculation was based on an 8-band model, included 
both two- and three-band terms in the two-photon am- 
plitude, but did not include terms with the so band as 
an intermediate stated More recently, we included the 
so band as an intermediate state, but only for two-band 
termsA4 The 2BTs in the 8-band model result differ from 
the 2BTs in the three-band model result of Sheik-Bahae 
by material independent factors. 

2. Spin Current 

The spin current PBA result is presented elsewhere** 
Here we summarize our earlier result in a new notation. 
For the electron spin current, [ic\ = MC2 = MC3 — 0, and 



m /m CiWl \ 3 / 2 m /to c , Wi \5/2 E P l-Z, 

HNl;e = D ; (l + Z c )-D 



m c \ to I ' m c \ to I 3E g 1 + xm c ,hh/mhh,ih 

m /'m C M l \ 3 / 2 {7 \ fm c , Ul \ b / 2 E P 1 



m c V m / \3 / m c V m / 



3E g 1 - xm c ,ih/mhh,ih 



m /"m Cthh \ 3 / 2 on m ( m c,hh\ b/2 e p l~Z, 

fJ-N2-e=D — — '■ — (l + Z c )-ZD 



ml m c \ m J 3E g 1 + xm Cl hh/™>hh,lh 

3 1 2 / 7 \ m ( m c ih \ 5 / 2 Ep 1 



D m rm^y/' U-lzA- Z C D^- 
Tfir- V m / V 3 / rrir V m J 



ml \ 3 / m c \ m I 3E g 1 - xm Ci i h /m h h,ih 



m fm c ,hh\ 5/2 e p l-Z, 

MJV3;e — —ZD 



m c V to / 3E g 1 + xm c j lh /m hh .i h 
m fm c , h \^/ 2 E P 1 



2 (1 - Z c ) D- 



m c \mJ 3E g 1 - xm c ^ h /m h h,ih 



(Bla) 



(Bib) 



(Blc) 



where x = (2Hu) - E g ) / (hui), m n } m = m n 1 - m" 1 , D is given in Ref0, and Z c = | E f£^ E ) ■ In ljBla|) and (|Blb|) 
{fJ>Ni;e and HN2-e), the first term is from the hh-c transition, the second term is from the hh-lh-c transition, the third 
term is from the Ih-c transition, and the fourth term is from the Ih-hh-c transition. In l|Blcfl for fJ,N^ e , the first term 
is from the hh-lh-c transition, and the second term is from the Ih-hh-c transition. Note that two-band terms make 
no contribution to /J.N3;e- 
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3. Spin 

To calculate optical effects due to the interference of allowed one-photon transitions and allowed-allowed two-photon 
transitions, we approximate the spin and velocity matrix elements and the energy denominator by their values at the 
r point, and approximate the energy bands in the (^-function as spherical and parabolic, neglecting the small k-linear 
term Ck and the small fc 3 spin-splitting. We used this method previously for two-photon spin injection. Since bands 
are degenerate at the T point, the lowest-order approximation to the matrix elements still depends on the direction 
k^oo However, by averaging the microscopic expression over physical systems rotated by each point group operation 
[which is equivalent to averaging over each term in Eq. I37al or Eq. I37b| , one can make the calculation using T-point 
states with pseudo-angular momentum quantized along z. The integral over k becomes a straightforward integral 
over the density of states in this approximation. 

The T-point basis states are given in (@}. However, all but the T§ c states are not eigenstates at the T point due 
to spin-orbit coupling between upper conduction and valence bands parameterized by A - . Using eigenstates to first 
order in A _ S we find 



o.=-^((^) 3/2 + (^n 



where 



Z+ = 



zi = - 



ClB = 

\J EpEpi 
A-Ep 



7T \ V TO / \ m / 



(2Hw) 
3/2\ ^J2hw - E, 



(B2a) 
(B2b) 



E' -Hw 



± 



A' - hw 



z'l 



3 

A-Ep, 



An 



K 



1 



1 



A' 



E' 



± 



Ao + fku 
1 



E' n +A E' - hw 



± 



E'o 



A' - hw 



i o 



hw 



In Z±, the first term is from intermediate sc states and the second term is from intermediate Ic and he states. In Z±, 
the first term is from intermediate so states, the second term is from intermediate sc states, and the third term is 
from intermediate Ic and he states. In Z±, the first term is from intermediate so states, and the second term is from 
intermediate he and Ic states. The term Z± can be neglected for typical semiconductors. Note that (Qa + 2(ib) has 
contributions only from intermediate so and sc states. This only includes transitions from initial hh and Ih states; 
transitions from initial so states, which contribute when 2hw > E g + A , have been neglected. 

4. Population 

We derive an expression for population control using the same method used above for spin control. To first order 
in A~. 



-e 3 2 

37T h 



V m J V to / 



3/2' 



v / 2hu> - Eg 
{2hwf 



^(X 1 +X 2 +X 3 ), 



(B3) 



where 



Xo 



-E r 



E' -hw E' + A' - hw 



2(E + A )- 1 (E'+A*,)- 1 f 2(E>+A Q )- 1 + (E>+A> 



E' Q -hw E' + A' - hw 



Aq + hw 



Xo = — 



A- E P , 



1 



1 



3 E' + A' \E' Q -hw E + A , -hw / 
Note that (— e 3 ) is positive. For typical semiconductors, X 3 can be neglected and 

A" 



Xo 



Xi 2 (A + hw) V Ep 



E F 



(B4) 
(B5) 
(B6) 
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In X2, the most important term is the last, which comes from the interference of {hh, lh}-so-c two-photon transitions 
and {hh, lh}-c one-photon transitions. 

The expression (IB3|) only includes the allowed-allowed transitions from the hh and Ih bands. At photon energies 
for which 2h\o > E g + Aq, one should add the contribution due to the transition so-uc-c. 

Because of <|21ll . (IF33I is also an analytical expression for Imx (2)a&c (-2u) ;lo,uj). Jha and Wynne have also used 
k-independent velocity matrix elements and spherical, parabolic bands to derive an expression for ^( 2 ) abc (— 2oj; oj, oj), 
but they did not include the interband spin-orbit coupling term A - ,* 4 - Taking the imaginary part of their Eq. 4.4 for 
hui < Eg < 2huj, and correcting a factor of 7r error, reproduces the lmx^ abc (— 2u>; u, (J) one would find from l|B3|) 
but with X2 — X3 — 0. Also, they make the approximation h\o rs E g /2 in the term X\. 

To get a PBA expression for the population control ratio requires PBA expressions for one- and two-photon 
absorption. We take the same approach used to derive (|B3(1 . but for simplicity, we take A~ = in the following. In 
the PBA, at photon energies 2tvu) < E g + Aq, one-photon absorption is 



'2mE F 



"W 3tt h 2 



V TO / V TO / 



y/2fuv - Eg 



(2huj) 2 



(B7) 



tv _ 
Hi 

In a material of cubic symmetry, the two-photon absorption tensor has three independent components 

£" 2 " bb , and £?2)° b , which are alternately parameterized by the set |C ( a 2 a ) aa ,cr, ( 5| (see Sec. US- The allowed-forb idden 
two-photon absorption in the isotropic Kane model, neglecting three- and four-band terms, is 



S(2) =?(2) 



2 2 / V fn V 6 6 



ki 



(B8) 



where 



f(2 



teV2e 4 E P (2hu-E g )* 



157T 



{2huf 



Note the additional symmetry, £?, 



aaa 
(2) 



2£2°P b + ffg) m this isotropic model. The allowed-allowed two-photon 



absorption, neglecting A' / (E' — E g + huj), has &§? a = & 2 ) 



abab 
(2) 



/-at 

4(i 



aabb _ o and 

- 2to E P ,E Q 



E P (E' -E g + hw) z 



which agrees with Arifzhanov and Ivchenko^S^ Thus, at photon energies for which allowed-allowed transitions domi- 
nate two-photon absorption, 



R 



to 



abba 
(2) 



(B9) 



whereas when allowed-forbidden transitions dominate two-photon absorption, 



R = 2hiu< 



EnEpi 



Ep (2hw - Eg 
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